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Abstract A brief physical-mathematical analysis of the electric stimulation of myocardium 
by current pulses delivered by pacing electrodes is done, introducing two hypothetical 
continua, electrically coupled, and partially superposed in space. From this analysis some 
desirable properties of the electric current field produced by the active electrodes are 
established. A theorem due to Prim is used as a guide to select possible shapes to be given 
to the active electrode heads to produce electric current fields with the established 
geometric properties. 
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(A) Introduction 


In the framework of artificial cardiac pacing, when the target region is close enough to the 
surface of the pacing electrode, there is evidence about a strong influence of the shape and 
size of electrode’s head on both, acute and chronic pacing threshold, and on the sensing 
performance of the electrode, ((G) 11, (G) 12, (G) 13, (G) 16, G (19)). 


' English translation of the 1984 paper “Teoria del electrodo de campo constante para la 
cardioestimulaci6n eléctrica crénica” with additional bibliography and comments. The 
original bibliography is cited as (G) followed by the corresponding reference number. The 
additional bibliography is cited as (H) followed by the corresponding reference number. 
Comments are given as footnotes. 
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The influence of the size and shape of the head of the electrode seems to be due mainly to 
the geometry of the electric current field generated by the electrode in the volume 
conductor 


formed by the biological tissues, and to the position of the target region (when pacing) or of 
the signal source (when sensing) in this volume conductor. 

All the procedures that aim to the electric stimulation of the heart can be understood using 
the following electrophysiological principle ((G) 7, (G) 8, (G) 14, G (16)): 


For a successful stimulation a critical quantity of electric charge must be transported 
through the excitable membranes of a suitable region (Lindeman’s liminal region) of 
myocardium depolarizing these membranes over a threshold membrane voltage, so that an 
action potential can emerge and propagate. 


To produce the above-mentioned charge transference, a transient electric current field 
(pulse) is produced by an active electrode located adjacent to the heart tissue that must be 
stimulated. To pace the heart a suitably spaced sequence of pulses must be applied. 

Now, let us consider the head of an intra-cavity electrode touching the right ventricular 
endocardium, at the end of a pacemaker lead, as shown in Figure 1. 


Figure 1: An active electrode is located at the tip of a catheter that is introduced (following a 
venous path) in the right cardiac cavity. 


Once the electrode is implanted, the presence of the foreign materials produces an 
inflammatory reaction followed by the formation of a capsule of fibrous tissue that 
surrounds the electrode’s head. The local properties of the volume conductor change. 
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As consequence, the so-called physiological increment of the stimulation threshold is 
produced ((G) 5, (G) 16, G (19)). Sometimes this ends chronic pacing due to an exit block: 
the pulse amplitude is too small to overcome the increased threshold. The capsule is 
sketched in the following Figure 2. 
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Figure 2: A sketch of the electrode lead with the electrode head separated from the excitable 
myocardium by fibrotic tissues. 


Besides, the relative movement between the electrode and the tissues, due to the contraction 
of the heart muscle, modifies the distance between the surface of the electrode head and the 
excitable myocardium. This distance variation is greater the stiffer is the catheter, and could 
produce an effect of mechanical hysteresis on the pacing thresholds. This mechanical effect 
could combine with an electric hysteresis effect (called Wedensky effect, see (G) 3). 

In the interface between the material of the active electrode head and the tissues, an 
electrochemical double layer appears. It is modified when an electric current pulse is 
applied by the artificial pacemaker. If the pacemaker produces a pulse of controlled 
voltage, the progressive polarization of the double layer during the pulse increases the 
impedance at the interface. As consequence the amplitude of the stimulating electric current 
pulse decreases. 


So, a system with strongly interacting subsystems is formed by the electrode and the 
tissues, whose electric, electrochemical, mechanical and biological properties must be taken 
into account in pacemaker electrode design. From a point of view that stresses a functional 
description of the signal flow through the system (like the one proposed in G (21)), it is 
convenient to distinguish four subsystems, partially superposed in space: ” 


> A more detailed description and discussion of these subsystems is given in H (8). 
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(a) A passive electric subsystem formed by four moduli connected in series: the 
electrode lead and head of the active electrode, the interface with biological tissues, 
the capsule of fibrous tissue that surrounds the electrode’s head, and the rest of the 
tissues ending in the indifferent electrode. Here the input is a current pulse that goes 
through each module between the electrodes. 


(b) An active electric subsystem formed by three moduli connected in series: the target 
myocardium near the active electrode, the rest of the excitable membranes of the 
cardiac muscle and the natural pacemaker. An electric signal output from the 
passive electric subsystem constitutes the input to the active electric subsystem. 
There are two outputs: an ECG signal and a physicochemical signal. 


(c) An active mechanical subsystem formed by the cardiac muscle receives the 
physicochemical signal from the active electric subsystem and contracts, producing 
a kinetic output. 


(d) A passive mechanical subsystem surrounding the active electrode’s head receives as 
input the above mentioned kinetic output from the active mechanical subsystem, 
and produces a kinematic output. This signal inputs in the passive electric 
subsystem, modifying the thickness of the capsule of fibrous tissue that surrounds 
the electrode’s head and varying the distance between the active electrode and the 
target myocardium. 


To simplify the mathematical modelling of this electrode-myocardium system, two facts 
must be considered: 


(1) There is a difference in order of magnitude in the electric resistance between of cell 
membranes from one side and of extracellular and intracellular compartments from 
the other side. 


(2) It is possible to distinguish four widely separated time scales in artificial cardiac 
pacing: 
(a) A scale of microseconds associated with the transients in tissues electrolytes 
reordering and inductive phenomena in the electrode lead. 
(b) A scale of milliseconds associated with the duration of pacing pulses and with 
the transport and capacitive phenomena that appear at the interface between the 
electrode and the biological tissues (electrochemical double layer). 
(c) A scale of tenths of seconds, that corresponds to the phases of cardiac cycle and 
to the mechanical phenomena in the electrode-myocardium system. 
(d) A scale of hours associated with the development of the tissues reaction to a 
foreign body (the electrode). 


To connect variables that belong to this time scales with variables of interest in a given 
reference scale, the variables in the slower scales can be considered as fixed parameters and 
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the variables in the faster scales can be considered relaxed to equilibrium in the scale of 
reference. > So, between a pacing pulse and the next one the distance between the 
electrodes and the target tissue could be modified due to muscle contraction. However, 
during the interval of application of a given pulse, in the scale (b), this distance can be 
considered as fixed, because all the modifications in the capsule of fibrous tissue that 
surrounds the electrode’s head occur in the time scales (c) and (d). 


(B) The two media hypothesis and the coupling between the two continua 


It is advisable to split the part of the electric subsystems belonging to the excitable 
myocardium in two superposed conductors of electricity: 

One conductor (let us call it the red continuum) is formed by the excitable membranes of 
the active electric subsystem. 

The other conductor (let us call it the blue continuum) is formed by the intracellular and the 
interstitial spaces of the passive electric subsystem combined in a volume conductor and 
electrically connected through the cellular membranes. Figure 3 shows an abstract 
representation of these continua. 


Figure 3: The two idealized continua. 


According to this approach, each point in the excitable myocardium belongs to both 
continua. + 


3 A more detailed discussion of these simplifying principles can be found in (H) 8, (H) 9 and 
(H) 10. 


* This was developed by Tung and other authors as the bidomain model of electrocardiology. 
(See H (6) and H (21)). 
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To obtain an operational (constructive) procedure of homogenization suitable to this case, 
the method of volume averaging already employed in the study of transport processes in 
porous media (G (18)) could be applied°. Meanwhile, a postulation approach allows write a 
pair of coupled equation for the red and blue continua in time tand space? . 


Considered as the abovementioned red continuum, the excitable membranes of the heart 
could be represented by a space filling three-dimensional nonlinear cable (a simplified 
isotropic and homogeneous version appears in G (8)). 

The blue continuum is anisotropic, and as shown in G (2), the anisotropy is different 
(unequal anisotropy) in the intracellular space relative to the interstitial space. This must be 
taken into account in the mathematical model. 


If V,, (t,7) is the field of membrane voltage (potential inside minus potential outside), a 
state variable of the red continuum, and if R,,R,,R, represent certain electric resistances per 


unit length in the direction of suitable chosen coordinate axis, assuming homogeneity 
(albeit unequal anisotropy) and other convenient assumptions, it seems possible to write the 
following approximate equation to connect the current density per unit volume across the 
excitable membranes /,, with an exogenous current density per unit volume /,. and with a 


density per unit volume of reactive currents associated with the non-uniformity of the 
polarization of the excitable membranes due to the active electrode: © 


2 2 
jE ee ean eR rec eal (1) 
R. Ox R, oy R, 02 


Currents are taken as positive when they flow from the intracellular space towards the 
extracellular (interstitial) space. 


> A foundation of bidomain model using the theory of averaged fields is given in H (9). A 
different foundation of the homogeneization procedure, using perturbation theory, is given in 
H (5). 

° Introducing the fraction y,, of membrane area per unit volume, the density of membrane 


current per unit area J, can be related with the density Sj (t,7) of electric current in the 


=Ve (C eV V,,)+ Ve (5 e J). As usual, V represents the 


tissue’s volume conductor by 7,, - J 


m 


gradient operator ande represents the scalar product. The tensor field C is symmetric and 
positive definite. It is related with a tensor generalization of the known scalar cable space 


constant. The symmetric tensor field D depends on the electrical properties of the blue 
continuum. The term V e(DeJ ) is equivalent to Rattay’s activating function of a nerve or 
skeletal muscle fibre. See H (9), H (10), H (18), H (19), H (20) and H (21). 
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We further assume that the blue continuum can be modelled as an anisotropic and Ohmic 
volume conductor. Applying the quasi-stationary equations of electromagnetics, it should 


be expected the following factorization: 
T(t?) =i(t)- F,() 2) 


; 0 ’ 
The membrane current is the sum of a displacement current C,, ae (here C,,, is an 


electric ae of excitable membranes per unit volume) with an _ ionic 


V 


m? 


current /,,,, W), where the vector W completes the set of state variables of the red 


continuum and embraces activation, recovery and other even slower variables related with 
the ionic channels. 

To simplify the analysis, it will be assumed that the activation variables are relaxed to 
equilibrium with membrane voltage, and the recovery variables will be considered as 


frozen. Then the ionic current /,,,, (V.) behaves as shown in Figure 4, being positive for 
membrane voltages between JV, (the rest membrane voltage, a first zero of the ionic current) 
and V,,,(the uniform membrane threshold, a second zero of the ionic current), and negative 
between V,,,and the excited state voltage V,(this last zero of the ionic current is a 


consequence of the freezing of recovery variables). 


Figure 4: Sketch of the ionic current per unit volume of the red continua (excitable 
membranes) as a function of membrane voltage for activation variables relaxed to 
equilibrium and recovery variables frozen. 


Considering the above mentioned two terms in the membrane current, equation (1) can be 
re-casted as follows: 
Ly oO L: : . 
Z oe ion V, )+ ° Z : : g y +1.,(t,7) (3) 
Zz 


m 
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(C) The treshold condition at membrane level 


After the end of a stimulating cathodic pulse the excitable membranes of a region located 
near the electrode are polarized due to the external current field. The space-time evolution 
of membrane voltages between stimulating pulses will be described by: 
é 1 @& L530 1 @& 
es 


Cp VV ey eee IY) 
mV re ax. m R, ay? m R, az? m (4) 


Ot m m 
According with Rushton, in the case of a nerve fibre, or according to Fozzard and 
Schomberg, in the case of a Purkinje fibre, in order to produce an action potential a liminal 
length of fibre must be depolarized above the uniform membrane threshold. A detailed 
study and justification can be found in G (8). 

Ventricular myocardium behaves as an electric syncytium. Lindeman’s experimental results 
about thresholds of stimulation of ventricular myocardium by external electrodes (G (14)) 
can be explained assuming that an action potential is produced if the membranes of a 
suitable three dimensional region are depolarized above the uniform membrane 


threshold V,,,. In principle the analysis given in Chapter 12 of G (8) could be extended 


x 


from 

a one dimensional nonlinear cable to a three dimensional one, and a derivation of the size 
and shapes of liminal regions obtained. ’ 

An heuristic approach to the liminal region could be this. Let us consider the more 


depolarized point with voltage Vj, and develop the ionic current /,,, (V,, ) in a 


neighborhood of V7, to the first order inVjy —V7y, . 


Applying a similar procedure as, for example, the one explained by Kyrala (G (10)) to get 
an interpretation of the Laplacian of a scalar field in a reference point in space, approximate 


: Vu-Vm a? 
the spatial second derivatives of voltage = >/m by — 2 pate “BY 
Ox ie Oy 
Vie —Vin a? Vin -V; 
~ and Dz. git m , all of them in the most depolarized point. The 


2 
ly 
lengths /,, / y and /, measure the distance between the most depolarized point and the 


points, located on the x, ¥, andz axes, where the membrane voltage is at its uniform 
threshold. They characterize the size and shape of an influence region of the electrode on 
the excitable tissue. 

Equation (4) is the approximated at the most depolarized point by: 

O 0 1 1 1 

Cees ee: en | Aa 5 
Ve pe le ah Ree <7] Var Yan) ©) 


ton (Vy, ) 


at aR, 


7 A full analytical approach to the characterization of liminal regions can be found in 
references H (9), H (10), H (14), H (18), H (19), and H (20). 
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: + : + 
2 2 
Ral, Bol ROE 


is greater (smaller) than 7 the 


Ov, 


m 


; wel 0 
From equation (5), if ee ion (Vin) 


voltage Vy, should increase (decrease). We could suspect that in the first case excitation 


will succeed and in the other it will fail. If the size of the influence region is too small, the 
excitation should fail. But all this lacks enough rigor. 


In any case we can guess that when the depolarization is more localized, its value in the 
most depolarized point should be greater, and the always negative term 
1: 1s, 16° Le 
-—5V,,+—: =, +—-= zl, Telated with the reactive currents there would be 
R, Ox R, Oy R, 0z 
greater in absolute value. As consequence, we could guess that other circumstances equal, 
the threshold cathodic pulse amplitude will be greater also. 
If this is so, it seems that to lower current threshold it would be advisable to depolarize 
as uniformly as possible the target region in excitable myocardium. An more uniform 
depolarization could be reached with a constant field electrode. To assess what could be 
attainable in practice, let us consider now the geometry of the current fields produced near 
the head of an active cathode. 


(D) Geometry of the exogenous electric current field in tissue’s volume conductor 


Let us consider a metallic cathode of any shape immersed in a volume conductor. 

Under quasi-steady conditions the field of electric current density J is solenoidal 
(divergence free): VeJ=0 (6) 
If ¢ is the tangent unitary vector to a field line [” in the point P then, with J = lV | being 


the local field magnitude (norm): 


=J-t (7) 
From (6) and (7) it follows: V i. 


J 
eJ=feVJ+J-Veti=0 (8) 


The gross point ® represents the scalar product of two vectors, V e( ) represents the 


” dJ . Sia Sete ; 
divergence operator and t e VJ = ads is the directional derivative of the current density 
Ss 


magnitude along the field line of intrinsic abscissa S (arc length) (G (1) and G (4)). 


Let us further assume that the volume conductor is isotropic. Then the unit vector t is at 
the same time tangent to a field line at a point P and a normal unitary vector to the 
equipotential surface crossed by the field line at the same point P (see Figures 5 and 6). 
There is a relation between the divergence of the field of normal unit vectors and the mean 
curvatures H of the equipotential surfaces at the same point (G (1), G (4) and G (9)). With 
a suitable sign convention for the mean curvature, it can be written thus: 


V ef(P)=-2-H(P) (9) 
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From equations (8) and (9) it is possible to show that current densities magnitudes J (P ) 
and J (0) at any two points Pand Q over the same field line I (see Figure 5) are 


related by a functional of the mean curvatures A r (also known as Germain’s curvature) of 
the points of the equipotential surfaces crossed by said field line between these two points: 


AB) ex ee (10) 


Figure 5: A planar sketch of field lines and Figure 6: A three-dimensional sketch of 
equipotential surfaces for a convex electrode. field lines and equipotential surfaces. 


The point O that appears in the interior of the cathode material in figure 5 is the electric 
centre of the electrode. The electric centre is such that taking it as origin of the position 


[(Jeii)-F-as=6 a 


OB 


vectors7 , the following equation is verified: 


In this formula the normal component J efi of the electric current density at a point in the 
surface of the electrode head is multiplied by the position vector of this point and the result 
is integrated on the entire electrode surface. 

From now on it is assumed that the conductivity of the electrode material (for example, 
Platinum) is at least an order of magnitude greater than the conductivity of the tissues. 

A good approximation is to consider that the conducting electrode surface is an 
equipotential surface and that the field lines are orthogonal to the electrode surface when 
the volume conductor behaves isotropically. 


10 
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Figure 7 shows a sketch of a concave electrode and its field of current lines, for an 
unbounded, homogeneous, and isotropic volume conductor. Some asymptotes are sketched 
also. 


¥% e Pectts c 


© centre 


Figure 7: The field lines [" are asymptotic to straight lines g that pass through the electric 
centre O of the electrode. The field line begins in a point S on the electrode surface that is 


in general different from the point S g Where the line 8 crosses the surface of the electrode. 


When the curvilinear abscissa S grows, the magnitude of the electric current density 
( ) i(t) 
JP) approaches the result corresponding to a point electrode: 
Pp Pp g Pp 4-2-(0S, ss) 


Let us introduce the excess mean curvature of an equipotential surface relative to a 
spherical surface constructed using the asymptote 2 : 


_ 1 


Bas ae (12) 


From Equation (10) and (12) we derive a formula for the magnitude of the electric current 
density in any point Pas a functional of the excess mean curvature integrated on a half 
field line beginning at P and extending to infinitum: 
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i(t) e 
Nae lis cae =e: BaG)e: (13) 


The generalization of formula (12) to the case of an anisotropic volume conductor will be 
considered elsewhere, applying some results obtained by Morris Kline (G (9)) to a problem 
of geometrical optics. However, the formula as it is by now should be enough for the 
purposes of this paper, centred in a qualitative discussion. 


(E) The possibility of a constant field electrode: the concept of an electrode’s critical 
region 


There is a theorem of classical vector field theory (due to Prim, see G (15)) that allows us 
to pose the type of electrode that is needed to make more uniform the polarizing effect of 
the cathode (the constant field electrode). 
According to it, for regular enough complex lamellar fields (fields that are equal to the 
product of a gradient field by a scalar field) the following three statements are equivalent: 

1- The surfaces orthogonal to the field are minimal surfaces (the mean curvature 

vanish at each point). 

2- The field is divergence free. 

3- The field magnitude is constant along each field line. 
A infinite Platinum half space electrode with a planar interface separating it from a 
complementary infinite half space filled with an ohmic, homogeneous and isotropic volume 
conductor produces a constant field everywhere, but is not attainable in practice. 
Nevertheless, a concave electrode like the one already sketched in Figure 7, shows several 
very interesting properties. 


Figure 8 sketches the behaviour of the magnitude of the electric current density J (P ) 


along a field line located near a cathode axis of symmetry, given as a 
function of the intrinsic abscissa s with the origin in the interface between the electrode and 
the tissues. 
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Figure 8: A sketch of the modification of the magnitude of the electric current density along a 
current line. 


12 


II International Congress of Biomathematics 


Buenos Aires, Argentina September 14, 1984 


The magnitude of the current density remains nearly constant along a segment of field line 


where the mean curvatures change their sign (the mean curvature in S, is exactly zero, 
being negative before and positive after it, according to the sign convention chosen here). 


An interval (s, AY) Jean be constructed such that inside it the mean curvatures remain in 
absolute magnitude less than a pre-established small number and the distance between the 
magnitude of the current density and its maximum value J (s.) remains less than a 


corresponding small number. 

Thus, it seems that a critical region will be produced, located at some distance from the 
cathode surface, where the field lines will be relatively parallel, and the magnitude of the 
current will be higher than in its neighbourhood. 


The assessment of the feasibility of designing, constructing and implanting an electrode 
with these advantages (“jumping” over the capsule of non-excitable tissues, lowering the 
electrochemical polarization of the double layer in the interface between the electrode and 
the tissues, and generating a critical region with suitable shape and size to lower the 
stimulation threshold counteracting the effect of the physiological increase) deserves a 
deeper study, both theoretical and experimental, with animal models and clinical trials. 


Merely as suggestions to begin with, Figure 9 shows two tentative electrode shapes. 
It is possible to anticipate that the one that appears in the left should show too much 
electrochemical polarization in the high curvature edges that surround the concavity. 
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Figure 9: Examples of possible electrode shapes. 
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Introduction The shape and size of electrode’s head have a strong influence both on acute 
and chronic pacing threshold and sensing performance when target’s region is close enough 
to electrode’s surface. This is particularly relevant to myocardial pacing, but besides that, it 
may be of interest for pacing of nervous system and skeletal and smooth muscle. 


Materials and Methods (a) Fields generated by seven model electrodes, of the same size 
and material composition but different shapes, were studied from an experimental point of 
view using an electrolytic tank. The equipotential surfaces surrounding each electrode were 
measured in steady state, for hemispherical, conventional cylindrical and five shapes of 
concave electrodes, in order to make a comparison of the resultant fields. (b) After that, a 
shape of concave electrode (named by us “clothoidal’”’) was selected and a real size platinum 
electrode was made with it. Then, it was implanted together with a conventional cylindrical 
platinum electrode, in the myocardium of ten dogs. The strength-duration curve for each 
electrode was measured in each dog, once a day, during several weeks. Rheobases, 
chronaxies and limit threshold charges were determined from experimental data. 


Results and Conclusions (a) A region located at a significant distance from electrode’s 
surface is formed, for special concave electrode (clothoidal), where current density is higher 
and current lines are relatively parallel (critical region of the electrode). (b) At implant, 
the strength-duration curve for the conventional electrode was under the strength-duration 
curve for the clothoidal electrode. But after several days the situation suffered an inversion, 
mainly because conventional rheobases began to grow significantly meanwhile clothoidal 
rheobases didn’t grew so much, maintained its acute values or even decreased. 


Key words: concave electrodes, clothoidal electrodes, critical region, stimulation 
thresholds, cardiac pacing, electrolytic tank experiments, scale models, rheobas, chronaxie, 
limit threshold charge. 
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EXTENDED SUMMARY 
(A) Introduction 


The shape and size of electrode’s head have a strong influence both on acute and 
chronic pacing threshold and on the sensing performance of the electrode, when the target 
region is close enough to the surface of the electrode (Irnich, 1980; Laplantif, 1983; Sudrez- 
Antola, 1983; Fiandra el alt, Chapters 7 and 8, 1985 a; Fiandra et alt, 1985 b; Aubert and 
Hector, 1985). This is particularly relevant to myocardial pacing, but besides that, it may be 
of interest for pacing of nervous system and skeletal and smooth muscle. 


The influence of the size and shape of the head of the electrode is due to the 
geometry of the electric current field generated by the electrode in the volume conductor 
formed by the biological tissues, and to the position of the target region (when pacing) or of 
the signal source (when sensing) in this volume conductor. 


In concave electrodes a certain critical region, where current density is higher and 
field lines are relatively parallel, is produced at a certain distance from electrode’s surface 
(Fiandra and Sudrez-Antola, 1987). The present experimental work affords some evidence 
that suggests an improvement of the pacing performance in a particular case of concave 
electrode, named clothoidal electrode, relative to other cardiac pacing electrodes. 


(B) Materials and Methods 


Fields of electric potential and current generated by seven giant model electrodes, of 
the same size (0.03 m diameter and 0.03 m height) and material composition (aluminium) 
but different shapes, were studied using an electrolytic tank (1 m diameter and 0.8 m 
height). 

The electrode’s shapes, all with cylindrical symmetry (in fact, flat meridian 
symmetry), are shown in vertical cross section in Figure 1: a hemispherical electrode, a 
plane electrode, two concave hemispherical electrodes (with arrows 0.004 m and 0.0085 
m), two clothoidal electrodes (with arrows 0.0137 m and 0.0170 m) and a hemispherical 
concave electrode with rounded borders. 

Figure 2 shows the measurement system: the tank (T) filled by an watery electrolytic 
solution (CV, homogeneous and isotropic volume conductor), the working electrode (R) 
mounted on an isolating cylinder (S) centred in the tank’s axis (e), the counter electrode 
(made of an aluminium foil of very large area), stick to the inner side of the vertical wall of 
the tank, a measuring electrode (E) (a slender pipette mounted in a mechanical systems that 
allows both horizontal and vertical displacements, and filled with a KCl solution 
surrounding a thin wire of Ag-AgCl), a voltage source and a voltmeter. A controlled (from 
0 to 9V) voltage source (F) was connected between the working and the indifferent 
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electrode as shown in Figure 2. The local voltage between the measuring electrode and 
earth was determined with a Hewlett-Packard multimeter (V). 


Figure 1: Vertical cross sections of the electrodes Figure 2: Electrical measurement system. 


A specific shape of clothoidal-concave electrode was chosen after analysing the 
results of the experiments done with the above mentioned electrolytic tank. 

Figure 3 shows in cross section the chosen shape for the clothoidal electrode used in 
the threshold experiments done in dogs. From this template, a special cutting tool was made 
in Barcelona, Spain, to be applied in a jeweler's lathe to obtain the desider concavity in one 
of the circular faces of a platinum cylinder. 


Figure 3: Template of a clothoidal electrode 


Real size clothoidal porous platinum electrodes were hand crafted by Dr. O. 
Fiandra (M. D.) in his jeweler's lathe. 

The clothoidal electrodes were implanted, together with a conventional cylindrical 
porous platinum electrode of the same size (2 mm diameter), in the myocardium of ten 
dogs. The strength-duration curves for each electrode were measured in each dog, once a 
day, during several weeks. Rheobases, chronaxies and limit threshold charges were 
determined from experimental data. Only some threshold values are reported here. The 
detailed experiments as well as a complete and in depth analysis of the obtained results will 
be presented elsewhere. 

During all threshold measurements, intracavity electrocardiograms were registered 
from both electrodes in order to compare their sensing capabilities. 
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Figures 4 to 6 shows one of the real size clothoidal electrodes implanted in a dog’s 
heart. 


Figure 4 Figure 5 Figure 6 


(C) Results 


(C.1) The equipotential curves sketched in Figures 6 toll were determined in steady 
state. We didn’t find significant differences in the equipotential surfaces determined with 
60 Hz AC relative to the same surfaces determined using DC. The equipotential curves of 
the hemispherical electrode are not shown because they correspond very closely to hemi- 
spheres at all distances from electrode’s surface, with the only exception of the region 
below the horizontal plane that passes through the electrode’s centre. This region has no 
interest, neither for pacing nor sensing. 


Beat: TSS 85) a= e a= 
eld] 


Figure 6 Figure 7 
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The electric potential difference between adjacent curves is the same in all cases, so that 
the reciprocal of the distance between curves gives a measure of the magnitude of the 
electric field. Due to the fact that the volume conductor behaves as homogeneous and 
isotropic, this reciprocal distance also can be taken as a measure of the local electric current 
density. 

As can be seen exemplified in Figures 7 to 11, in the concave electrodes case, a 
particular region centred on the axis of symmetry and located at a certain distance from 
electrode’s surface is formed. For all electrodes the electric potential decreases 
monotonically along the axis of symmetry. But for concave electrodes the potential shows 
an inflection point where the magnitude of the electric field (and as consequence, the 
magnitude of the electric current density) is a maximum. In a neighbourhood of this point 
the magnitude of the electric current density is higher and the current lines near the axis are 
relatively parallel. 

This region can be named the critical region of the electrode. This critical region is 
both, farther away (from electrode’s surface) and bigger, for clothoidal electrodes in 
comparison with other concave electrodes of the same size and material. 
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(C.2) Now we show the rheobase, the chronaxie and the limit threshold charge calculated 
from the measured strength-duration curves in two dogs (D-4 and D-5) both for the 
clothoidal (a) and the conventional (8) porous platinum electrodes. The data correspond to 
the implantation of the electrodes (time 0), and then to the first (time 1), the fourth (time 4), 
and the eight (time 8) day after the implant. 


D-4 
Time 0 Time 1 Time 4 Time 8 
Rheobase a 0.20 0.71 0.71 0.72 
(mA) B 0.49 0.68 1.12 1.13 
Chronaxie o 0.54 0.36 0.21 0.38 
(ms) B 0.35 0.15 0.29 0.30 
Limit Threshold «a 0.11 0.25 0.22 0.27 
Charge B 0.17 0.10 0.32 0.34 
D-5 
Time 0 Time 1 Time 4 Time 8 
Rheobase a 0.64 1.22 0.84 0.81 
(mA) B 0.39 1.05 1.55 2.15 
Chronaxie a 0.27 0.21 0.17 0.49 
(ms) B 0.21 0.24 0.15 0.31 
Limit Threshold a 0.17 0.26 0.14 0.40 
Charge B 0.08 0.25 0.23 0.67 


The rheobases behave according to the physiological increment of the threshold of 
excitability of myocardium (Fiandra and others, 1985 a, Chapter 7; Fiandra, Sudrez-Antola 
and others, 1981). 

In D-4 and D-5, as well as in the other dogs, the rheobases of the clothoidal 
electrodes increased less than the rheobases of the conventional electrodes. The rheobase of 
the clothoidal electrode in D-5 at Time 0 is well above the rheobase of the conventional 
electrode, but at Time 8 the clothoidal rheobase is well below the conventional one. The 
limit threshold charge shows a similar behaviour. This was a common finding in our 
experiments, as will be reported elsewhere. In D-4 and D-5, relative to their values at 


implant, chronaxies of both, clothoidal and conventional electrodes, first decrease and then 
increase. 


(D) Conclusions 
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(a) For a concave electrode a region located at a certain distance from electrode’s 
surface is formed, where current density is higher and current lines are relatively 
parallel (critical region of the electrode). 


(b) This critical region is both, farther away (from electrode’s surface) and bigger, for 
clothoidal electrodes in comparison with other electrodes of the same size and 
material. 


(c) At implant, the strength-duration curve for the conventional electrode was in general 
under the strength-duration curve for the clothoidal electrode. But after several days 
the situation suffered an inversion, mainly because conventional rheobases began to 
grow significantly meanwhile clothoidal rheobases didn’t grew so much, maintained 
its acute values or even decreased. 


(d) In general, the clothoidal electrode showed a better sensing performance relative to 
the conventional electrode implanted in the same dog heart. 
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Introduction The geometry of the electrical current field generated by a given electrode in 
the volume conductor formed by biological tissues and the properties and position of the 
target region in excitable tissues, both determine the performance of the electrode as a 
stimulating and as a sensing device. Both from experimental facts and from theoretical 
reasons, the geometry of the current field of certain concave electrodes deserves special 
attention. 


Materials and Methods For a field with a symmetry axis the distribution of current 
sources at electrode’s surface is substituted by multipolar coefficients assigned to the 
electric centre of the electrode. Using Levenberg-Marquardt algorithm, the multipolar 
coefficients can be determined from suitable experimental results. Besides this, the current 
densities at any two points over the same field line are related by a functional of the mean 
curvatures of the points of the equipotential surfaces crossed by said field line between 
these two points. Then, Prim’s theorem for solenoidal and almost gradient fields can be 
applied to the analysis of the critical region of the electrode, experimentally found as 
described in the first part of this work. 


Results and Conclusions At points not too far away from the symmetry axis, the field can 
be well represented by three or four terms of the multipolar expansion, even in the case of 
concave electrodes. The electric potential, considered as a function of distance on the 
symmetry axis, shows an inflection point for concave electrodes and none for convex 
electrodes. This inflection point is the centre of the critical region of the electrode. The 
critical region of the electrode can be defined as a certain neighbourhood of the surface 
formed by the points where the equipotential surfaces have zero mean curvature. Its 
position and extension is the related with the multipolar coefficients of the field. 


Key words: concave electrodes, clothoidal electrode, multipolar expansion of a potential 
field, solenoidal fields, almost gradient fields, critical regions, mean curvature, Prim’s 
theorem. 
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EXTENDED SUMMARY 
(A) Introduction 


The geometry of the electrical current field generated by a given electrode in the 
volume conductor formed by biological tissues and the properties and position of the target 
region in excitable tissues, both determine the performance of the electrode as a stimulating 
and as a sensing device. Both from experimental facts and from theoretical reasons, the 
geometry of the current field of certain concave electrodes deserves special attention. 


Figure 1. Sketch of the equipotential curves, Figure 2. Sketch of the equipotential curves, 
clothoidal concave electrode hemispherical concave electrode 


(B) Materials and Methods 


(B.1) Figure 3 shows a sketch of an electrode with an axis of rotational symmetry. 


Figure 3. Sketch of an electrode’s head with a set of polar coordinates rand @ on each 
symmetry plane through the axis Z . 
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The distribution of current sources at electrode’s surface can be substituted by a set 
of multipolar coefficients assigned to the electric centre of the electrode. Then for an 
unbounded, isotropic and homogeneous volume conductor, the electric potential V(r, @) can 
be developed in terms of Legendre’s polynomials and reciprocal powers of the polar 
distance r , beingc, the 2” multipolar coefficient (Kellog, 1929): 


n=+00 


vir,a=> = - P,(cos 6) (1) 


n=0 


When the volume conductor is bounded but the symmetry of revolution is retained, it 
is necessary to add a constant term d, and terms involving increasing powers of I . 


The origin of coordinates Q (see Figure 3), when the electrode injects a net current 
in the volume conductor, can always be chosen to coincide with the so called electric 
centre, so that the coefficient of the dipolar termc, =0 . Taking this into account, we found 
that it is possible to approximate the equipotential surfaces measured in the upper region of 
the electrolytic tank, with focus in a neighbourhood of the z axis, by a formula that 
includes only a constant, a monopolar, a dipolar and an octupolar terms: 


V(r,0) = d, ++ 2- P,(cos0)+ =: B(cos 6) (2) 
r r r 


The unknown multipolar parameters as well as the position of the electric centre in 
the symmetry axis can be estimated minimizing the following chi-squared variable, where 
V(7,,9,) represents a measured voltage value at point (7,,0,) in a vertical plane through 


the symmetry axis: 


2 
N 

F aeese)= S [v.89 —d, — © _ 2. B(cos 0,)- ; P,(cos 6, ) (3) 
k=l KON, NK 


In this case k = 1,2,..., N represents the points where the measurements were done. 
The distance r, between the electric centre of the electrode and the point where the voltage 


is measured is given by: he = ly. +P, +2-1 +P, COS QD, (4) 


In formula (4): (1) , 1s the distance between the point 0 , of intersection of the symmetry 
axis z with the surface of the electrode’s head, as shown in Figure 3, and the point P, where 


the voltage is determined; (2) @, is the angle between the segment 0, P, and the z axis; 


and 7, is the distance between (Q and 0 f (see Figure 3). It and can be called electric 
radius of the electrode. While p,and g,can be measured directly, (7,,0,) must be 
determined after the electric radius has been estimated, so in practice we have to work with 
the operational coordinates (p,,@, ). 

To estimate the multipolar coefficients and the electric radius 7, it is possible to 


apply the Levenberg-Marquardt algorithm (Marquardt, 1963; Press and others, 1992) in two 
steps. First a gross estimation is obtained working with the voltage values on the electrode 
axis only. In this case we put @=0 andr=Z+1%, so the following formula is used to 
estimate the unknown parameters 7, dy), Cy.C55C3: 
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V(z) = d, spp he ae 
eth (z+m) (etn) 

Then, this gross estimation is used as seed to obtain a more accurate estimation employing 

all the measured values of voltage, according to formulae (3) and (4). 

If A,is the area of the conductive electrode’s head, then to begin the nonlinear iterations 


A, 


rcan be equalled to = A seed for parameters d,, cy c, and c, can be estimated 
‘Te 


sequentially (as done in compartment analysis) from the numerical values of V(z) when the 


(5) 


distance on the axis of symmetry, from the surface of the electrode, grows. 


(B. 2) Let us consider now a so called almost gradient field, like the electric current 


density J in an ohmic but perhaps heterogeneous volume conductor of variable (from 
point to point) scalar conductivity G : 


J=-G-VV (6) 
If the vectorial field is solenoidal, by definition the divergence of the field vanishes: 
VeJ=0 (7) 


According to Prim’s theorem (Prim, 1948; Sudrez-Antola, 1984) for solenoidal and 
almost gradient fields, two of the following three conditions imply the third one: 
1- The orthogonal surfaces to the field have zero mean curvatures. 
2-The field is solenoidal. 
3- The magnitude of the field along a field line is constant. 

Equation (7) can be applied to almost stationary fields, that is, in conditions of slow 
enough time variation to be able to work with static equations. This is the case commonly 
encountered while pacing the heart or during functional electric stimulation of nerve and 
muscle fibres (Fiandra and others, 1985 a, Chapter 8; Reilly and others, 1992). 

However, the volume conductor formed by biological tissues is often anisotropic. 
This is particularly important in the case of ventricular myocardium (Zipes and Jalife, 1990, 
Chapter 20). Now the relation between the current density vector field and the gradient of 
electric potential field is given by the equation 


J=-G-VV (11) 


The electrical conductivity G is a positive definite symmetrical tensor. The current density 
field is not exactly almost-gradient as previously assumed. The bidomain model of 
ventricular myocardium, that is now used to analyse threshold and action potential 
propagation, takes into account the so called unequal anisotropy (different anisotropy ratios 
in the intacellular domain in comparison with the extracellular domain) and has important 
practical consequences (Sudrez-Antola, 1994 a) 


(C) Results 


(C.1) Figures 4 and 5 show the adjusting of formula (2) on the axis of symmetry of both, 
the plane and the clothoidal concave electrode (arrow 1,37cm), respectively, minimizing the 
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corresponding chi-squared by Levenberg-Marquardt algorithm. Formula (5) was employed 
with c, fixed at zero. 
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Figure 4 Figure 5 
We found: 


Plane electrode: %=1.05cm d, =1.66V Cy =2.13Vxcem c, =-0.92V xcm’* 
Clothoidal concave: 7 =0.87cm d,=0.84V cc, =1.01Vxcm_  c, =-4.76 Vxam’* 
Since the multipolar coefficients c, are proportional to the injected current, a more 


meaningful set of parameters is given by Z, = oy 
Co 
For the clothoidal electrode with arrow 1,70cm, we found, including now the octupolar 
coefficient 8, (and with a v* =1,88x10~): 
B, = 2 =-0.466 cm? B; = =0.089 om’ 
Co sa) 
For a clothoidal electrode with arrow 1,37cm, we found adjusting parameters with the 
complete set of measured voltages, including also the octupolar coefficient /;, : 
By = 2 =-0.65 cm’? B; = =0.25 em? 
Co 0 
A complete report of the estimation of parameters will be given elsewhere. 


(C.2) Now, consider the mean curvatures of the equipotential surfaces near the axis of 
symmetry of a concave electrode in an isotropic and ohmic volume conductor. Near the 
concave surface of the electrode, mean curvatures are positive. But if we move along the 
axis away from the surface of the electrode, the mean curvature decreases until it vanishes 
at a certain point on the symmetry axis. In this point the curve that represents the voltage as 
a function of the abscissa Z has an inflection point. This can be seen in figures 4 and 5. 

In the case of a clothoidal electrode this inflection point is fairly distant from the 
electrode. In the case of a plane electrode, the inflection point doesn’t exist. 

Now, it is possible to apply Prim’s theorem to a concave electrode taking into due 
account the continuity properties of the current field. We see that in the neighbourhood of 
the surface with zero mean curvature, the field lines will be almost parallel, the magnitude 
of the current density will be nearly constant along the field line and the magnitude of the 
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current density will be higher than its magnitude in points of the same field line away from 
the surface with zero mean curvature. 

As consequence, a critical region is formed near the axis of a concave electrode, 
where the field lines are nearly parallel and more densely distributed. 


(D) Conclusions 


At points not too far away from the symmetry axis, the field can be well represented 
by three or four terms of the multipolar expansion, even in the case of concave electrodes. 

The electric potential, considered as a function of distance on the symmetry axis, 
shows an inflection point for concave electrodes and none for convex electrodes. 

This inflection point is the centre of the critical region of the electrode. 

The critical region of the electrode can be defined as a certain neighbourhood of the surface 
formed by the points where the equipotential surfaces have zero mean curvature. Its 
position and extension is the related with the multipolar coefficients of the field. 
Figure 6 shows a qualitative picture of the aforementioned critical region. In the last article 
of this series (Sudérez-Antola and Artucio, 1994) the position, size and shape of the critical 
region will be studied quantitatively. 


Figure 6 ‘Figure 7 — Figure 8 


Figure 7 and 8 shows a template that was used to construct an optimum concave electrode 
and a photograph of the clothoidal electrode hand crafted by Dr. O. Fiandra using his 
jewelry’s lather. 
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CONCAVE ELECTRODES UI: COMPUTER 
ASSISTED DESIGN 


Roberto Suarez Antola and Gabriel Artucio 
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Introduction The design of concave electrodes able to generate a critical region (C.R.) 
where current density is higher and field lines are relatively parallel, located at a certain 
distance from electrode’s surface may be of interest both for electrical simulation of living 
tissues and for sensing tissue’s own electrical activity (1). We present here a draft of a new 
design method that allows a sequential determination of size and shape of electrode’s head, 
such that the electrode will be able to generate a C.R. of pre-established position and shape 
in given volume conductor, assuming an axis of symmetry for the desired field. 


Materials and Methods The connection between geometric parameter of electrode’s head 
and its C.R. is established trough the multipolar coefficients of the electric potential field 
(2). In a first stage we obtain an approximation to the monopole, quadrupole and octupole 
coefficients as functions of geometric parameter of C.R. (This relation is obtained from 
Whittaker’s classical solution of Laplace’s equation). In a second stage we simulate the 
potential field by means of a finite analytical approximation and make comparison of 
obtained C.R. with the desired one. If necessary we add some multipoles and proceed to a 
new digital simulation, and so on until we reach the desired degree of accuracy. 


Results and Conclusions Then, we identify the point on the axis of symmetry of the field, 
where the potential has its extreme value before going into monotone variation towards 
zero. The equipotential surface that crosses the axis in this point can be adopted as the best 
model for size and shape of electrode’s head. The proposed method doesn’t involve the use 
of ab-initio numerical algorithms for the solution of the exterior boundary value problem of 
potential theory. 


Key words: concave electrodes, clothoidal electrode, size and shape of the critical region, 
Whittaker’s solution of Laplace equation, multipole expansions. 
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EXTENDED SUMMARY 
(A) Introduction 


The shape and size of electrode’s head have a strong influence both on acute and 
chronic pacing threshold and on the sensing performance of the electrode, when the target 
region is close enough to the surface of the electrode (Sudrez-Antola, Griego y Fiandra, 
1994). 


The influence of the size and shape of the head of the electrode is due to the 
geometry of the electric current field generated by the electrode in the volume conductor 
formed by the biological tissues, and to the position of the target region (when pacing) or of 
the signal source (when sensing) in this volume conductor. 


When the volume conductor can be considered as approximately homogeneous, 
isotropic and unbounded, the axisymmetric field of electric potential generated by a given 
symmetric electrode can be considered as produced by a certain set of multipolar moments. 
If the multipolar are located at the so-called electric centre of the current distribution the 
dipolar moment is always zero. At points not too far away from the symmetry axis, the field 
can be well represented by three or four terms of the aforementioned multipolar expansion: 
monopole, quadrupole, octupole and so on (Suarez-Antola, 1994). 


This is so even in the case of concave electrodes. In concave electrodes a certain 
critical region, where current density is higher and field lines are relatively parallel, is 
produced at a certain distance from electrode’s surface (Suarez-Antola, Griego y Fiandra, 
1994; Suarez-Antola, 1994). We present here a sequential method to design a concave 
electrode’s head such that it will be able to generate a critical region of pre-established 
position and shape in a given homogeneous and isotropic volume conductor, assuming an 
axis of symmetry for the desired potential field. 


(B) Materials and Methods 


Let us begin with Whittaker’s solution (Copson, 1975) for Laplace’s equation in 
cylindrical coordinates, written for a field with the axis z as its symmetry axis: 


2a 
Vir.)=5— | fe+i-1 0086) dO [1] 
0 


Here r is the distance between the considered point and the axis. When r=0 (that is, 
on the axis z), V(0,z)=f(z), so that given the electric potential over the axis of symmetry its 
value at any other point can be found by formula [1]. 
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Then we develop V(r,z) in powers of r up to order four at least (the coefficients of 
odd powers vanishes), first by approximating f(z+ircos®) by a Taylor polynomial in ircos® 
and then by calculating the error from Whittaker’s formula. 


The next stage is to introduce the multipolar expansion along Zz axis: 


co 


I : : : : ge 
(z= ee B ~_ where Lis the total electric current, G is the electric conductivity of the 
n=0 Z 


tissues and the 4 are the multipolar coefficients, being A, = 1. If the origin z=0 is at the 
electric centere, then the dipolar coefficient #,=0 and after the monopolar term we have to 


consider the quadrupole 3, the octupole B4 and so on. 


In order to generate a critical region as was mentioned in the introduction, 83 must 
be negative. Then, at a certain distance from z=0, along the z axis, f(z) (for I positive) 


reaches a maximum value for a certain z=zy and then tends monotonically towards zero 
2 


while z grows. For z=Z,, f(z) shows an inflection point (that is ae f(z.) =0). The point 
ce 
(r=0, z=Z.) 18 going to be the centre of the intended critical region. We then consider 


E.=-ZVr0 and E=-LV(0,9). When r=0, Er is always zero. E, =|E. (0, z) is 
Zz r 


maximum for z = Zc. 


le 


Cc 


Then we consider the region whose points (r,z) are such that 1—- ¢< <1l+e and 


IE, 


— <é fora given ¢ suitable chosen (for example e=0.1). 


Cc 


Using a polynomial approximation for V(r,z) with a finite multipolar expansion for 
f(z), we can determine the size and shape of the critical region as a function of the field 
parameters. Conversely, given the dimension of the intended critical region, it is possible to 
estimate the corresponding field parameters. The size and shape of the surface of the 
electrode that generates such a critical region can be determined constructing the equi- 
potential surface that passes through the point z=zm™ of the axis of the field. If the result thus 
obtained is not good enough, another multipolar coefficient can be taken into account and 
the whole procedure must be repeated again. 


(C) Results and Conclusions 
The equipotential surface that crosses the point of the axis after which the fields 


decreases monotonically can be adopted as the best model for the size and shape of 
electrode’s head. 


Physics in Medicine and Biology, vol. 39 a, 1994 


The proposed method doesn’t invoke the use of ab-initio numerical algorithms for 
the solution of the exterior boundary value problem of potential theory. 


Instead, it only invokes the solution of fairly simple and well conditioned nonlinear 
algebraic equations (Press and others, 1992). After a certain size and shape for the 
electrode’s head is selected, a detailed numerical study of the field can be done by more 
conventional means, or a sound physical scale model can be constructed and the potential 
field can be measured. 
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Sensing bioelectric signals with concave electrodes: 


I. Lead Fields and Sensitivity 
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Abstract 
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The shape and size of electrode’s head may have a strong influence on the sensing performance, both of intra- 
cavitary electrodes and intra-tissue electrodes, whenever the sensed region is close enough to electrode’s 
surface. The purpose of this paper is to discuss the performance of concave electrodes when used as sensing 
devices. The shape and size of a critical region of higher sensing sensitivity to the electric activity of biological 
tissues is studied using a mathematical model of a homogeneous and isotropic volume conductor. A formula is 
derived for the measured potential difference between a proximal concave electrode and a distal indifferent one, 


due to tissues’ electrical activity. 


Keywords: bioelectric sensors, bioelectric signals, concave electrodes. 


(A) Introduction: 


Sensors of bioelectric phenomena are usually 
known as electrodes. As every other biomedical 
sensor, they serve as interfaces between biological 
tissues and electronic equipment for diagnostic or 
research purposes. Both biologic and electronic 
parameters have influence on _ the sensing 
performance of the electrode, so that they must be 
considered during electrode’s design. Sensing 
electrodes are designed to fulfil with three 
conditions: (1) Provide an_ output — signal 
proportional to the measured bioelectric quantity. 
(2) To be insensitive to other physical quantities. 
(3) Minimize the unavoidable effects due to 
bioelectric activities of other origin and noises of 
several kinds. Intracavitary and _ intratissue 
electrodes are often used both for diagnosis and 
therapy, so that the two kinds of application must 
be considered during electrode’s design. 


This is the case of intracavitary electrodes used to 
take cardiac electrograms and to pace the heart. 
This is done with the same electrode located at the 
tip of a catheter that is introduced (following a 
venous path) in the right cardiac cavities. Once the 
electrode is implanted, the presence of foreign 
materials produces an inflammatory reaction 
followed by the formation of a capsule of fibrous 
tissue that surrounds the electrode’s head. 
Therefore, the local properties of the volume 
conductor change, so that the sensing and the 
stimulating properties of the electrode may change. 
This and many other problems must be considered 
in the process of design of heart electrodes, most of 
them with convex shapes but with geometric details 


that are introduced to improve their sensing and 
pacing performance [1]. 


During 1983 appeared the idea of pacing the heart 
with concave electrodes, instead of the then 
standard convex ones. Some prototypes of concave 
electrodes were developed in France [2] (parabolic- 
concave), in Uruguay [3] (clothoidal-concave), and 
later in other places. 


Now, it is well known that both the shape and size 
of electrode’s head have a strong influence on acute 
and chronic pacing threshold as well as on the 
sensing performance of the electrode if the target 
region (for stimulation or sensing) is close enough 
to electrode’s surface. Although this influence is 
particularly relevant in the case of myocardium 
pacing [4], it could be of interest both for pacing 
[5] and sensing the local electrical activity of the 
nervous system, skeletal muscle, and smooth 
muscle. The sensing properties of concave 
electrodes have not been much studied, and even 
then, only experimentally and in a global approach 
[6], without considering the distribution of local 
sensitivities as given by the lead field of the 
sensing electrodes. 


The lead field concept arises when two electrodes 
are placed in a volume conductor forming what is 
called a lead [7], [8]. 


The lead vector (7) describes how an 
electromotive source of vector density E',(t,7) 


(located at the point of position vector 7 and at the 
time instant ¢) influences the lead voltage V, 


> 


recorded between the electrodes, being u-E, a 
scalar product: 
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Vy, t= [arog t,r dv (1) 
B 
The integration is, in principle, extended to the 
whole  tissue’s ohmic volume conductor 
represented by the three-dimensional region B . 


For a given lead vector, the contribution of a given 
vector density of electromotive force is 
proportional to the cosine of the angle between 


them. If E s is orthogonal to uw it will not 


contribute to the measured voltage signal, and the 
contribution will be maximum in absolute value if 
the two vectors are parallel. The contribution of a 
given local electromotive force density will be 
proportional to the magnitude of the local lead 
vector (7 ). So, the local lead vector describes the 


sensitivity of the lead to a source located there. 
The distribution of the lead vectors in the volume 
conductor is the lead field, and this field gives the 
distribution of sensitivities of the system (the lead) 
formed by the volume conductor (biological 
tissues) and the electrodes. It follows from the 
reciprocity principle [8] that the lead field (7) 


is identical with the field of electric current 
density that is produced in the volume conductor 
when a unit current (called the reciprocal 
current) is fed to the lead from the outside. The 
directions of the lead field flow lines give the 
directions of the sensitivity, and their density is 
proportional to the local magnitude of the 
sensitivity. This density is higher near each 
electrode, so that ceteris paribus (other things 
remaining the same) the electric signals produced 
in the neighbourhood of the electrode receive a 
higher weight in the determination of the output 


signal V, t . For convex electrodes the sensitivity 
grows when the distance to the electrode decreases 


and reaches a maximum value at the interface 
between the electrode and the tissues. For example, 


in the case of a spherical electrode of radius 7,, 


immersed in an isotropic, homogeneous and ohmic 
volume conductor, the lead field is given by: 


me 
Je, G8 (2) 


2 
4nr 


ur,d,0 =( 


Here ris the distance from the centre of the 
electrode, and @,@are the angular spherical 
coordinates of the unit vector é, directed outwards. 


The magnitude of the lead field vectors is 
maximum on the surface of the electrode, where 


rreaches its minimum value7,. (Equation (2) 


follows from the relation 477° j vr =T between 
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the global current / injected through the electrode 
and the current density 7, taking into account that 


the magnitude of the local lead vector is 


ur = F and due to the assumed spherical 


symmetry is directed radially outwards). 


The purpose of this paper is to discuss some 
interesting properties of concave electrodes when 
used as sensing devices. This is done in the 
framework of the theorem of reciprocity for a 
heterogeneous and anisotropic conductive medium 
with a distribution of current sources 
(electromotive forces) in its interior and a pair of 
sensing electrodes on its boundary. One electrode 
is active, concave, and located very close to the 
region of interest, and the other indifferent, with a 
much bigger surface, and far from this region. 


(B) The lead field for concave electrodes 


Now we will show that for concave electrodes of a 
suitably selected shape there is a critical region, 
located at a certain distance from the surface of the 
electrode, where the density of lead field lines is 
higher than in any other region (including regions 
located nearer to electrode’s concave surface) and 
the lines are relatively parallel (Fig. 1) . 


Fig.1: Pictorial representation of a meridian cross section of a 
concave electrode. Sketch of the lead field lines and of some 
equipotential surfaces. The critical region appears at some 
distance adjacent to the axis of symmetry of the lead field. 


The shape and size of this critical region of higher 
sensing sensitivity to the electric activity of 
biological tissues may be studied using a 
mathematical model of a homogeneous and 
isotropic volume conductor, as functions of the size 
and shape of the surface of the electrode. 

So, we consider a symmetric electrode immersed in 
an unbounded ohmic volume conductor, with the 
indifferent electrode located far away. 
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In this case the electric potential field inside the 
volume conductor verifies Laplace’s equation. A 
convenient starting point in this case is Whittaker’s 
solution for Laplace’s equation [9], assuming 
cylindrical symmetry in the electric potential field 
V inside the volume conductor: 


20 
V p,z = ff z+jpcosé dé (3) 
one, 


Here p is the distance between the considered 
point and thez axis. When p=0 (that is, on 
the z axis) (0,z)= f(z), so that given the electric 


potential over the axis of symmetry its value at any 
other point can be found by formula (3). 


Then we develop V(p,z) in powers of p up to 
order four at least (the coefficients of odd powers 
vanish), first by approximating f/f z+ jpcosé 
by a Taylor polynomial in jpcos@ and then by 
calculating the corresponding terms and the error 
R(p,z) from Whittaker’s integral formula: 
V(p,z)=Zo(z)+ °Z2(z)+p*Z4(z)+R(p,z) (4) 
By definition: 
1 d?Z,(z 
Z(2)= fl), Za(2)=—4 0b 
4 dz 
ld “Zale) 

Zz 5 

4(z) 64 ast ( ) 


Furthermore, the remainder: R(p, z)= o(° ) 


The next stage is to introduce the multi- 
polar expansion of the electric potential along 


_1 24, 
ae 7G onl (6) 


Z axis: Liza i 


In this formula I is the total injected electric 

current, G is the electric conductivity of the tissues 

and the parameters B are the so called multi- 
n 


polar coefficients. 3, =1. In the case of a potential 


field with a nonzero mono-polar coefficient, the 
origin z =Qcan be suitably chosen at the electric 
centre of the distribution of field sources [10]. Then 
the dipolar coefficient £,=0 and after the mono- 


polar term we must consider the quadrupole ,, the 
octupole £, and so on. The following expansion for 
the potential V(0,z) is enough for our purposes: 


I \l 6. B 
zae—|-+34+4 7 
f 4nG|z z z* @) 
The local lead vector sensitivity is calculated from 
: OV oV 
its scalar components, u, = uy = for a 
a °° ap 
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unit global current / =1 , using the first multi-polar 
coefficients , and f, as the only parameters. 


yw | VO)_ 2? PO), 470] (94) 
i dz 4 dp 64. dz 


.|_e VS) , 30° a* Fle) 
u,® a ae + a as (8b) 


The coefficients 6, and f, can then be related 


with the size and shape of the sensing electrode’s 
head. This can be done either measuring the 
distribution of the electric potential adjacent to the 
symmetry axis and estimating the parameters from 
the experimental data [6] or making a digital 
simulation of the potential field produced by the 
electrode. To generate a critical region as was 


mentioned in the introduction (, must be negative 
(Fig. 2). 


f(z) 


Fig.2: Pictorial qualitative sketch of the voltage on the axis of 
symmetry of the sensing concave electrode. 


Then, at a certain distance from z=0, along the 
z axis, f(z), for/ positive, (in the lead field’s case 


equal to 1), reaches a maximum value for a certain 
z=z,,and then tends monotonically towards zero 


while z grows. The function f(z) has a negative 
slope for z greater thanz,,. When z is near enough 
to Z,, its second derivative is negative. However, 


for z positive and big enough, the function tends 
to zero with positive second derivative. Then, there 
is a certain z=z, such that f(z)shows an 


2 
inflection point: SS f(z.) =9 
Z 


The point 9 =0,z =z, is going to be the centre of 
the intended critical region. We then consider the 
scalar components of the lead vector, according to 
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formulae (8a) and (8b). When p =0, then the radial 
component wu, is always zero. The absolute value 
of the axial component u,,;.(z)=|w-(0,z)| takes its 


maximum for z =z,. 


Now we have: 


3 
He)~ se.) LEN e 2.) API. y 
i" (z.) 
Let us consider the length scale: /,” = ; si (9) 
fe.) 
dz* 


This scale characterizes the length along the z axis 
in which the component of the lead fieldu, , suffers 
a significant variation. 


Then, we define the dimensionless space 
coordinates: 


Z 


=> (10a) n= (10b) 


Cc Cc 


From equations (8a) and (8b), and from the 
definitions (10a) and (10b), it follows that the 
components of the lead field near the point: p =0 


z=z, can be expressed as: 


u 
Me pw (Sl) (11a) 
Uu. 2 
u & " 
pet ae oregon Lb 
u, ag 
3 
Here, by definition, s = vey HE 
Z 


It is equal to 1 in our case. Then we consider a 
region extending symmetrically around the z axis, 
and centred in z =Z,. whose points 9,Z are such 
that ge eres and bol < for a given ¢é 
Uc Uc 
suitable chosen (for example ¢ = 0.1). This region is 
a critical region of the leading field in the sense 
mentioned at the beginning of this section (B). The 
sensitivity lines are almost parallel and their 
density there is greater than in other regions of the 
lead field. This is so because the magnitudes of the 
local lead vectors are very near their possible 
maximum value u, in any region adjacent to the 


Z axis. 
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(C) The lead voltage: a general formula 


Concave electrodes show variable curvatures, in 
magnitude and in sign, from point to point on their 
surfaces. The concavity is surrounded by a convex 
ring. The current density is smaller in the points 
located on the concavity than in the points located 
in the convex portions [11]. As consequence, the 
degree of polarization of the interface is greater in 
the convex part than in the concave one, so that the 
voltage drop will be different at different points of 
the electrode’s surface. The influence of the voltage 
drop in the electrochemical double layer located at 
the interface between the concave electrode and the 
tissues, for a polarizing concave electrode, may be 
taken into account by means of the introduction of 
a family of virtual non-polarizing concave 
electrodes of smoothly varying shapes. In each 
instant, one electrode belonging to the family 
represents the local time behaviour of the 
polarizing electrode. The voltage drop at the 
indifferent electrode is assumed negligible due to 
its big area. The geometry of the equivalent volume 
conductor is identical with the real one, except in 
the vicinity of the already mentioned family of non- 
polarizing concave electrodes. 


Let us sketch the derivation of a new formula that 
generalizes formula (1) to this more complex 
situation. 


We represent by Gthe conductivity tensor of the 
tissues, and by ie -G-VV+G-E, the current 
density produced by the combined effect of a 
potential gradient and a local electromotive force. 


This is the current field produced in the volume 
conductor by the electrical activity of the tissues. 


Under quasi-stationary conditions, the density 

current field generated by the tissues verifies: 
V:-J=0 

The conductivity tensor is symmetric, so that 

a:G=G-a for any vector d. 

In the boundary of the volume conductor there 

aren’t electromotive forces and with the only 


exception of the points that belong to the surfaces 
of the electrodes, the boundary condition is 


J-n=0. 
On the electrode surfaces the quasi-stationary 


electric potentials can be taken as constants and 
equal to V,and to V, respectively, so that the lead 


voltage is given by: V, =V, -V, 
When the current is injected through the lead from 
outside, the current density will be represented 


by 7 =-G-Ve , being ¢ the induced electric 
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potential in the tissues due to the external current 
only. 


We also have Vi =( and the same boundary 
conditions for ¢ as given before for V . 


Now, we consider stationary alternating currents of 
frequency: w=27f 


The phasor amplitudes of the _ scalars 


potentials w,7 , V W,P and the phasor 


=> = 


amplitudes of the vector fields J CF Oe 


verify the following equations: 
Ve J=VWV-j—j Es (12) 


From vector analysis tools it follows: 

J VV -jdv= | V- Vi dv—{V Vij du 

B B B 

Considering that V oy =0 we obtain: 

fvViae=fV ij dv=V,1 (13 a) 

B OB 

Here [ = { jen dS is the phasor amplitude of 
dB, 

current introduced through the electrodes into the 


tissues volume conductor. 
Also, we have: 


vad ae oJ Pay V-J dv 


Considering that V-J =0 we obtain: 


fv: i dv=[¢ Rey dv= 4-6, I’ 
B OB 


(13 b) 


I' is the phasor amplitude of the global current 
produced by the electromotive forces associated 
with the electric activity of the biological tissues. 


From all these relations we derive 


[Ve-J dv = pe. T= 
B 
VI jE, dv 

i 


(14) 
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Then, if Z, 


xt 


w is the impedance due to the 
external measuring circuit between the electrodes, 


Vow 27. (15) 


If Ze « is the impedance due to the in-series 


combination of the ohmic resistance of the 
tissues R,.with the lumped impedance Z, w of 


the polarizing interface between the electrode and 
the tissues, then: 


bw—-gw -Z,wlw (16) 
ae ee : 
Introducing u =F Jj. from equations (14), (15), 


and (16) it follows, after several calculations: 


: 1 + > > > 
Vi w =— (fi WF +E, Wr di) 
(= B 
Loy W 


(17) 
This formula is a nontrivial generalization of 
equation (1). The phasor component V, w of the 


voltage sensed by the lead that corresponds to the 
frequency @ , is expressed as a functional of the 
distribution of phasor vector electromotive 


forces, E, w,r , in the tissue’s volume conductor, 


and the phasor lead’s vector field @ w,7 


produced in the volume conductor by an alternating 
unit current injected from outside through the 


sensing electrodes. The lead field W w,7 may 


depend on the frequency due to the different 
degrees of polarization of the points on the surface 
of a concave electrode if the mean curvatures differ 
between these points. 


(The impedance Z, w may be given by the 
formula derived in [13]: 
R, 


Wy 


Bie 


Here mis a positive number, less than 1, related 
with the fractal dimension of the interface, R, is a 


D.C. resistance of the polarized interface, and @, is 


a characteristic angular frequency that separates the 
constant phase angle behaviour of the interface at 
high frequencies from its resistive faradic 
behaviour at low frequencies). 
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If the impedance of the external circuit is at least an 
order of magnitude greater than the combined 
impedance of the tissues and interfaces, (17) 


simplifies into: V, w = fii wr -E, w,r dv 
B 


When the polarization of the interface does not 
modify significatively the lead field, 


Uw, eur 
Then, returning to the time domain we obtain 
formula (1) given in the introduction. 


(D) Discussion and conclusions 


(a) The problem of determining the geometric 
parameters of an electrode’s head such that it 
produces a critical region with suitable shape and 
size can be approximately solved by an analytical 
approach. Let us observe first that the electric 
conductivity of the electrode’s material is at least 
several orders of magnitude greater than the 
conductivity of the tissues and the body fluids. As 
consequence, the surface of the electrode may be 
safely assumed to be a surface of equal potential of 
the electric potential field produced by the 
electrodes when the reciprocal current is injected to 
the lead from outside. In principle any surface 
V(p,z)=constant could be a candidate to a surface 


of a suitable electrode. However, this cannot be 
exactly this way if the surface of equal potential is 
unbounded. But if our interest is focussed into the 
abovementioned critical region located near the 
axis of symmetry, it is possible to construct a 
suitable approximation to the shape and size of the 
electrode to produce the desired field within 
established tolerance bounds. This will be studied 
in the second part of the present work. 
Nevertheless, the procedure can be briefly sketched 
for an unbounded, isotropic and homogeneous 
volume conductor. 


Let us consider Whittaker’s formula again. We 
introduce a function z=g(p) such _ that: 
Qn 


1 ; 
Veep =— ff g p +jpcosd dd=V, 
20% 
(18) 
The function z=g(p) represents the trace of the 
surface of constant potential,V,, in a plane that 


contains the axis of symmetry of the field. In 
principle, given V, and assuming enough 


regularity in the functions f(z) and g(p) we can 
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calculate all the derivatives of g(p) ing =0, from 
(18), as functions of the parameters #, and £, 

This can be done by successive derivation 
operations under the integral sign, relatives top, 
and putting each result equal to zero. 


An approximation of order four may thus be 
obtained: 


d°g 0 F 1d‘g0 , 
2 dz’ 24 dz" 


Z=S p 2+ 
(19) 


For a concave electrode of depth d and radius Ry, 
we may write: 
zag p xm t— op —-So' 
R, R. 


0 

(20) 
The geometric parameters d and R, are determined 
from the derivatives of the function g()on the 


axis of symmetry, and these derivatives are then 
related with the multi-polar coefficients of the field. 
The parameter z,must be chosen near but 


somewhat greater than the parameter z,, where the 
function f(z)takes its maximum value, in order 


that the current density may have the same sign 
everywhere in the surface of the electrode. 


Then z.-d will give the distance, along the 


symmetry axis, from the centre of the critical 
region to the electrode. See Fig. 3. 


AZ 


Fig.3: Qualitative sketch of the curve corresponding to formula 
(20). The depth d and the radius Rp of the concave electrode 


appear, as well as z ,, and the centre of the critical region z,.. 


(b) The lead field adjacent to a polarizing concave 
electrode was studied from a mathematical 
viewpoint, and the resultant distribution of sensing 
sensitivities was given through approximate 
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analytical formulae for the scalar components of 
the lead field (8a) and (8b). 


(c) The critical region for sensing with concave 
electrodes was defined and approximate analytical 
formulae for the scalar components of the lead field 
were obtained in (1 la) and (1 1b). 


(d) A generalization of formula (1), that relates the 
lead voltage V(t) with the lead field a()and the 


distribution of electromotive forces E s £7 inthe 


volume conductor, was given, for alternating 
current signals, in formula (17). The amplitudes 
ii(@,7) of the oscillating lead field may differ 


significantly from u(*) only near the concave 


electrode surface. Far from this surface the effects 
due to the curvatures vanish and the potential field 
resembles the field of a convex spherical electrode. 


(e) The lead field was studied using a model of an 
isotropic and homogeneous volume conductor. In 
the cases of many tissues, and in the case of 
myocardium in particular, the volume conductor is 
anisotropic and heterogeneous, so the study should 
be extended to embrace this more general case. 


(f) Every concave electrode has curved convex 
borders. The magnitude of the lead field (and as 
consequence the sensitivity to bio-electric signals) 
is higher there than in the critical region. 
Nevertheless, the fibrous capsule of connective 
tissue that surrounds the chronic pacing electrodes, 
in most cases extends over this region and beyond. 
So, this potentially higher sensitivity has no chance 
to manifest itself after a few weeks from the 
implantation. 
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Abstract 


As shown in the first part of this work, the shape and size of electrode’s head may have a strong influence on the 
sensing performance, both of intra-cavitary electrodes and intra-tissue electrodes, whenever the sensed region is 
close enough to electrode’s surface. The purpose of this second paper is threefold: (a) Review some 
experimental results obtained for the lead fields and the corresponding sensitivities of a set of metal electrodes 
of different shapes and suitably enlarged size to use an electrolytic tank and comparing these results with the 
mathematical model for concave electrodes immersed in homogeneous and isotropic volume conductors 
developed in the previous paper. (b) Derive in the time domain a formula for the measured potential difference 
between a polarisable proximal concave electrode and a distal indifferent one, due to tissues’ electrical activity. 
(c) Give some guidelines to be applied in the conceptual design of concave electrode heads intended to sense 
bioelectric signals generated in electromotive sources imbedded in heterogeneous and anisotropic tissues. A 
virtual space (alibi), homogeneous and isotropic, is introduced in correspondence with physical space. The 
design problem is easier to pose and analyze in this auxiliary space. Then the results may be transformed back 


to the physical space. 


Keywords: bioelectric sensors, bioelectric signals, concave electrodes, tissue’s anisotropy, and heterogeneity. 


Introduction: 


Sensors used to detect and measure bioelectric 
phenomena are known as electrodes. They serve as 
interfaces between biological tissues and electronic 
equipment, for diagnostic, therapeutic or research 
purposes. Both biologic and electronic parameters have 
influence on the sensing performance of the electrode, 
so that they must be taken into account during 
electrode’s design. 


Sensing electrodes are designed to fulfil with three 
conditions: 


(1) Provide an output signal proportional to the 
measured bioelectric quantity. 


(2) To be insensitive to other physical quantities. (3) 
Minimize the unavoidable effects due to bioelectric 
activities of other origin and noises of several kinds. 


As was shown in the first paper of this series [1], both 
the shape and size of electrode’s head may have a 
strong influence on the sensing performance of an 
electrode. This occurs if the target region for sensing 
(that is the region where the distribution of 
electromotive forces of interest is located), is close 
enough to the surface of the sensor. 

Although the shape of most electrode heads is convex, 
several years ago a generation of concave electrode 
heads was introduced, mainly but not exclusively, with 


the intention to improve electrode’s performance during 
cardiac pacing [1], [2], [3]. To characterize the sensing 
properties of concave electrodes, it is necessary to study 
first the distribution of local sensitivities in the tissue’s 
volume conductor as given by the lead field [4], [5] of 
the sensing electrodes. The (for the time being assumed 
stationary) lead vector u(F) describes how an 


electromotive source of vector density E (tr), 
varying with time ¢ and located at the point of position 
vector 7 influences the lead voltage V,(t) recorded 
between the electrodes: 


ViO= fa@)-EwAdv (a) 
B 


The integration is, in principle, extended to the whole 
tissue’s ohmic volume conductor represented by the 


three-dimensional regionB (u.E, is a scalar 


product). For a given lead vector, the contribution of a 
given vector density of electromotive force is 
proportional to the cosine of the angle between them. If 


E, is orthogonal to uw it will not contribute to the 


measured voltage signal, and the contribution will be 
maximum in absolute value if the two vectors are 
parallel. The contribution of a given local electromotive 
force density will be proportional to the magnitude of 


the local lead vector U (7 ). So, the stationary local lead 
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vector describes the sensitivity of the lead to a source 
located there. The distribution of the lead vectors in 
the volume conductor is the lead field, and this field 
gives the distribution of sensitivities of the system (the 
lead) formed by the volume conductor (biological 
tissues) and the electrodes. 


It follows from the reciprocity principle [4] that the 
lead field ii(7) is identical with the field of electric 


current density that is produced in the volume 
conductor when a unit current (called the reciprocal 
current) is fed to the lead from the outside. In our case 
one electrode will be active, concave, and located very 
close to the region of interest, and the other indifferent, 
with a much bigger surface, and far from this region. 


In reference [1], assuming cylindrical symmetry in the 
electric potential field inside the volume conductor, 
approximated analytical formulae for the components 


of the lead vector field ii(7) were obtained for any 


electrode’s head size and shape, with the only 
restriction of the abovementioned imposed cylindrical 
symmetry and the assumed homogeneity and isotropy 
of the volume conductor. This was done starting from 
Whittaker’s solution for Laplace’s equation [6], using a 
multipolar expansion [7] of the electric potential along 
the axis of symmetry of the field. 


It was shown that for concave electrodes of a suitably 
selected shape there is a critical region, of maximum 
sensitivity for sensing purposes, located at a certain 
distance from the surface of the electrode, where the 
density of lead field lines is higher than in any other 
region (including regions located nearer to electrode’s 
concavity) and the lines are relatively parallel (Fig. 1). 


Also, a generalization of formula (1) was derived for 
the measured potential difference between a proximal 
concave electrode and a distal indifferent one, due to 
tissues’ electrical activity [1]. This formula, valid for 
steady alternating currents of a given (but arbitrary) 
frequency, may be applied even if the volume 
conductor is heterogeneous and anisotropic. 


This second paper aims to: 


(a) Review, from the standpoint of the lead field and the 
corresponding sensing sensitivities, some experimental 
results obtained in an electrolytic tank with metal 
electrodes of suitably enlarged sizes and made of a 
material showing strong polarization of the metal- 
electrolyte interface. Compare the results of the 
measurements with the mathematical model for 
concave electrodes immersed in homogeneous and 
isotropic volume conductors developed in [1]. 


(b) Derive for the time domain a generalized formula 
that gives the measured potential difference between a 
polarizing proximal concave electrode and a distal 
indifferent one, due to tissues’ electrical activity, for 
transient’s currents and voltages. 


(c) Extend the design guidelines of concave electrode’s 
heads (for attaining critical regions of suitable locations 
and sizes), briefly suggested in [1], to embrace the 
selective sensing of bioelectric signals generated in 
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electromotive sources imbedded in heterogeneous and 
anisotropic tissues. 


Fig.1: Pictorial representation, taken from [1], of a meridian cross 
section of a concave electrode, with a sketch of the lead field lines 
and of some surfaces of constant potential. The cross section of the 
critical region appears (in gray) as fuzzy elliptical. 


Measurements of the lead field for concave 
electrodes 


The surfaces of constant potential surrounding each 
type of a set of seven models of electrode heads of the 
same size but different shapes were measured in an 
electrolytic tank of 100 cm diameter and 80 cm height, 
with an electrically isolating internal surface (painted) 
[8],[9]. 

The indifferent electrode was made with a symmetric 
band of an aluminium foil stick to the inner side of the 
vertical wall of the tank. 


The electrode heads (a hemispherical one, a planar one, 
and five different shapes of concave ones) were lathe 
machined on the top face of aluminium cylinders of 3 
cm diameter and 3 cm height. 


Each cylinder was mounted on an isolating (plastic) 
cylindrical support of 5 cm diameter and 5 cm height 
and located in the centre of the isolated bottom face of 
the tank. 


The top face and the upper part of the curved vertical 
face of each cylinder were in electrical contact with the 
electrolyte, and the rest of the curved face was 
electrically isolated. 


An Ag/AgCl measuring electrode was mounted in a 
vertical slender pipette, filled with a suitable 
electrolytic solution. Then the measured voltage may be 
assigned to the point where the tip of the pipette is 
located. 


A controlled voltage source (dc and ac) was used to 
produce the fields in the electrolytic solution filling the 
tank. The measurements were made using a digital 
multimeter of very high input impedance. Also, an 
oscilloscope was used to register transients during 
sensing experiments with short current pulses. Figure 2 
shows a sketch of the experimental arrangement. 
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Fig. 2: Draft of the experimental arrangement. R, electrode; S, 
support; E, pipette with Ag/AgCl electrode; D, mechanical dispositive 
to locate the pipette; T, tank wall with the inner surface isolated and 
with the indifferent electrode stick to it. 


In this paper we consider the measured lines of constant 
potential, in a plane through the axis of symmetry, for 
only two electrodes immersed in the electrolytic 
volume conductor of the tank: one representative of the 
concave electrodes (named clothoidal according to 
[3]), the other flat. 


The measurements of the lines of equal potential were 
done with d.c. and once a steady state was established 
in the system. See figures 3 and 4: 


Fig.3: Distribution of lines of constant electric potential for a 
concave electrode in steady state. Difference of potential between 
adjacent lines: 0.05 V 
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Fig.4: Distribution of lines of constant electric potential for a flat 
electrode in steady state. Difference of potential between adjacent 
lines: 0.10 V 

The corresponding orthogonal lines, giving a graphical 
representation of the lead field and its sensing 
sensitivities, could be approximately constructed from 
the lines of constant potential applying the classical 
drawing recipes of the old engineering field theory. 


Instead of doing that, we use a mathematical model to 
connect the experimental results with the theoretical 
approach developed in [1]. 


Let us remind the formula for the potential near the axis 
of symmetry (where the critical region of highest 
sensing sensitivity of the concave electrodes is located: 


V(p,2)=Z,(z) + Z,(z) + p'Z,(z) + R(p,z) 
(2) 
Here p is the distance between the considered point 


and the axis of symmetry, and Z is the coordinate on the 
axis. 


With V(0,z)= f(z): 


Pod? A (z) 
2a re re te 
64 dz 
The following expansion for the potential V(0,z) 
should be enough for our purposes: 


fe~a++ S454 


3 
4rG|z Z @) 


The remainder verifies: 


R(p,z) = 0(p") 

In the present case/is the global electric current 
injected by the voltage source in the steady state, and G 
is the electric conductivity of the electrolyte in the tank. 
The origin Zz =Ovwas suitably chosen at the electric 
centre of the distribution of field sources, so that the 
monopole coefficient 8,=0. In practice we only know 


that this centre must be on the axis of symmetry. So, let 
us choose the origin of coordinates in the point of 
intersection of the z axis with the upper face of the 
cylinder of the electrode, and measure distances d 
uppwards along the axis from this known point. Then 
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we putz=r7,+d, being7yy)a new parameter, the so 


called “electric radius” of the electrode. For a perfectly 
spherical (and non-polarizable) electrode the geometric 
radius and the electric radius are equal. 


After the mono-polar term, with 6)=1, we involve 
only the quadrupole and the octupole terms. 

The multi-polar coefficients 8, and £,, as well as the 
electric radius7), must be determined from the 
experimental results. 


The parameterc did not appear in the multipole 
expansion of the potential given in [1]. This expansion 
is valid when the volume conductor, besides being 
homogeneous and isotropic, may be considered infinite. 
The additional parameter must be introduced to take 
into account that now the volume conductor cannot be 
considered as unbounded. In general, the smaller the 
geometric dimensions of the electrolytic tank for a 
given electrode, the greater will be the number of 
additional parameters and additive terms necessary to 
quantify boundary effects on the spatial distribution of 
the electric potential. 


In our case a single additive parameter is enough 
because: the geometric dimensions of the aluminium 
cylinders (radius 1.5 cm) are an order of magnitude less 
than the smallest dimension of the volume conductor 
(tank radius 50 cm, tank height 80 cm); the upper face 
of the electrode is located 8 cm from the isolated 
bottom of the tank; the indifferent electrode covers a 
significant part of the inside of the vertical wall leaving 
a relatively small isolating band in the lower part of the 
wall, band that begins slightly below the level of the 
electrode head (Fig. 2); and, finally, the electric 
potential measurements were taken mainly near the 
electrode head. A thorough discussion of boundary 
effects on the distribution of the potential in the 
electrolytic tank may be found [9]. 


Let us consider N values of the potential measured on 
the axis: 


V(0,7r, +d,)=V,,n=1,2,...N 


Then, a set of values of7),f,,f,andcmay be 


obtained minimizing the nonlinear expression 


2g» Bos By) = SF (ly + dys By» Buse) —V,) 


(4a) 
ae +d,,8,,B3.C) © 
rea iee 8, b, 
¥ i on 
anGin+d, | “Ht+dy “Htay 
(4b) 


This was done both for the concave and the flat 
electrode using the powerful Levenberg-Marquardt 
algorithm [10]. 

The following results were obtained from 15 
measurements of the potential on the axis of symmetry 
for each electrode: 
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Clothoidal electrode: 

ro 1.00cm, By ~ —0.53 cm? B, ~0.01cm*, 
c~1.00cm'! 

Flat electrode: 

ro 1.05 cm, 8B, ~ —0.40 cm? B, = 0.00 cm}, 
c~1.66 cm! 


In both cases the final Chi-squared values were of the 
order of 10%. The adjustment is very good. In fact, 
leaving aside, for obvious reasons, the hemispherical 
electrode, it was also very good for the other four types 
of concave electrodes of different shapes. 


Three remarks are indicated to help to interpret the 
experimental results and the estimation of parameters 
presented above: 


(a)The conductive surface in contact with the 
electrolytic solution (real exposed area) was 
approximately the same for the seven types of 
electrodes studied with the electrolytic tank. 


(b) No difference was found between the surfaces of 
constant potential measured with a sustained a.c., 60 
Hz, relative to those measured with d.c. in the steady 
state, although the measured points were not very near 
to the electrode surface. 


(c) The polarization and other effects of the indifferent 
electrode can be neglected in this case. 

Now, let us turn to the problem of determine the field of 
local sensitivities from the experimental results in the 
electrolytic tank. As was shown in [1], the local lead 


vector (sensitivity) W (F ) may be calculated from its 


scalar components, axial wu, = _ and radial 
Iz 
oV ‘ : 
u,, =-—— for a unit global current / =1 , using the 
P 6p 
first multi-polar coefficients f, and £; as the only 
parameters. 
uw] Ue) 2 VP) 0° PFE)) Gay 
: dz 4 dF 64 dz 
2 3 44 
d 
ipo lee2 f(z) , 30 le) (50) 
2 dz’ 64 dz 


This presupposes that the location of the electric centre 
is known, which in practice, is not the case. However, 
using formula (3) with z=7, +d , we introduce a new 


unknown parameter (the electric radius) and the known 
distance (measurable) from the electrode head. 


The parameterc, necessary to achieve a _ good 
adjustment between formula (3) and the measured 
values of the potential along the z axis, disappears from 
the formulae (5) because the components of the lead 
field are functions of the first and higher derivatives of 


(3). 
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The lead voltage: a formula in the time 
domain for sensing transient signals with 
polarisable electrodes. 


Formula (1), reminded in the introduction, can be 
obtained under the following assumptions: 


(a) The module of the input impedance of the 
measurement chain is, for all the frequencies involved 
in the biological electromotive sources, at least an order 
of magnitude greater than the module of the in-series 
combination of the impedance of the electrode’s 
electrochemical double layer with the impedance of 
biological tissues. 


(b) The effect, on the lead field, of the transients and 
non-uniformities in the polarization of the double layer 
located at the interface between the surface of the 
electrode head and the electrolytic volume conductor, 
may be neglected. 


Note that when a direct current is established in an 
electrode-electrolytic system, (as is the case of the 
electrodes in the electrolytic tank studied in the 
precedent section), once the polarization of the double 
layer reaches a steady state, the lead vector field in each 
point of the volume conductor remains constant. 

If the conductive medium is homogeneous and 
isotropic, the electric potential field may still be 
described by formulae (2) and (3) (withz=% +d). 
The steady and non-uniform polarization of the double 
layer appears implicit in the values of the 
parameters 7), 8,, and f; . 


In principle, if these parameters are given as suitable 
functions of time, the effect of transients and non- 
uniformities in the polarization of the double layer, on 
the axial u, and radial u, components of the lead 


field, could be described. 

The non-uniform polarization of the double layer in 
non-spherical electrodes is related with the non-uniform 
distribution of the electric current density in the surface 
of the electrode, due to differences in the mean 
curvatures between points on the surface. 

This non-uniform distribution is suggested in Figure 5 
for a convex and non-polarisable electrode. 

There, a continuous variation of the mean curvature is 
suggested in a cross section of the electrode. The point 
O represents the electric centre of the electrode. 

The points S, P and Q are located on the same lead field 
line. 

The dotted line between points O and S is a virtual 
prolongation of the lead field line. 

It is possible to derive an analytical formula relating the 
magnitude of the lead field at each point of a field line 
with the curvatures of the congruence of orthogonal 
surfaces crossed by the line, by a suitable 
reinterpretation of results developed in [11]. 


p 
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Fig. 5: Sketch of the distribution of the lead field lines and the 
surfaces of constant potential near a convex electrode surface 
immersed in a homogeneous and isotropic volume conductor. 


This non-uniform distribution is suggested in Figure 5 
for a convex and non-polarisable electrode. There, a 
continuous variation of the mean curvature is suggested 
in a cross section of the electrode. The point O 
represents the electric centre of the electrode. The 
points S, P and Q are located on the same lead field 
line. The dotted line between points O and S is a virtual 
prolongation of the lead field line. It is possible to 
derive an analytical formula relating the magnitude of 
the lead field at each point of a field line with the 
curvatures of the congruence of orthogonal surfaces 
crossed by the line, by a suitable reinterpretation of 
results developed in [11]. 


When an electrode surface is in equilibrium with an 
electrolytic solution, a steady and uniform potential 
difference appears between the two phases. If a direct 
current is injected through a highly polarizable flat 
electrode with a sharp edge, the potential drop in the 
double layer is modified: the dynamic perturbation is 
greater where the curvature is greater. At the fist 
moment after the current is applied, the perturbation is 
everywhere negligible, and the lead field (assumed 
quasi-static) corresponds to a truly flat electrode 
uniformly polarized. But as time goes on, a distribution 
of different potential drops in different places on the 
surface of the electrode begins to develop, and the 
corresponding lead field changes as it were produced by 
a family of non-polarizable virtual electrodes that 
change progressively their shapes, until arriving to a 
steady state. Each non polarizable virtual electrode may 
be characterized by its electric radius and its multipolar 
coefficients. The whole process can be described with a 
suitable set of time dependent parameters. Figure 6 
suggests what happens once the steady state is arrived: 


Fig. 6: Flat electrode strongly polarized in steady state, together with 
the equivalent virtual electrode in steady state. 


The inner surface of the virtual electrode is practically 
coincident with the surface of the flat electrode in 
points where the curvature of this last one is zero (flat 
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upper face far from the edge) or relatively small 
(vertical face far from the edge). On the edge and near 
it, the virtual electrode is rounded and appears as 
smaller relative to the real flat one. 

When the volume conductor is heterogeneous and 
anisotropic, the approach through the electric radius and 
the multipole coefficients doesn’t apply. In [1] the 
generalization of formula (1) for the measured voltage 
signal in the time domain was done in the frequency 
domain. 

For steady alternating current signals of frequency @ , 
finite input impedances of the chains of measurement 


Ze.(@) and non homogeneous and non isotropic 


volume conductors, taking into account the combined 
impedance Z, (a) of biological tissues and_ that 


produced by the polarization of the interface between 
the electrode and the tissues, the following formula 
was derived: 


Lo 


tJ tlw?) E(w?) dv) 
(i+——) 3 
Z,, (Ww) 


ext 


(6) 


The component V, (w) of the voltage sensed by the lead 


that corresponds to the angular frequency @ , is 
expressed as a functional of the distribution of vector 


electromotive forces, E',(w,7), in the tissue’s volume 


conductor, and the lead’s vector field u(w,7) 


produced in the volume conductor by an alternating unit 
current injected from outside through the sensing 
electrodes and in steady state of oscillation. As 
mentioned above, the lead field may depend on the 
frequency due to the different degrees of polarization of 
the points on the surface of a concave electrode if the 
mean curvatures differ between these points. 


Formula (6) is useful mainly when the spectral content 
of the signals doesn’t change significantly with time. 
But in biomedical engineering there are many 
circumstances in which the spectral content of relevant 
signals changes significantly during the interval of 
measurement, and a mixed time-frequency description 
may be necessary [12]. 

So, to allow the application of the powerful modern 
tools of signal analysis, a derivation in the time domain 
is advisable. 


Let us begin with the general relation, equivalent to 
Equation (14) of [1], but written in the time domain: 


(O-pOlO= 


Se = = = 7 
VOU + f HG?)-EtF) dv Se 
B 
Equation (7) can be derived from a version of the 
reciprocity theorem formulated for the so-called 
quasistatic approximation to the electromagnetic 
field. 
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To simplify the derivation, let us suppose that the ohmic 


component of the external circuit impedance, R,,,, is 


dominant. Then, the voltage signal (lead voltage 


between the electrodes) V,(t) and the current 


I' produced by the bioelectric electromotive forces are 
related by: 


V, (t) x R.. I (t) (8) 
The injected current ys produced by an external applied 


potential difference between the electrodes Y, —Y , 


is assumed constant: the applied potential difference 
must be adjusted to comply with this restriction while 
the electrode’s electrochemical double layer polarizes. 


Then, if R, is the resistance of the volume conductor 


t 
of the tissues and [ Kw du is the voltage drop in 
0 


the double layer: 


Pr —Y, =|R, + | K(u)du I (9) 
0 


(In [13] may be found a short presentation of this kind 
of mathematical model of the electric behavior of 
polarizable electrodes in volume conductors). 

From (7), (8) and (9), ifR 


exziS at least an order of 


t 
magnitude greater than R, + { K(u) du, follows the 


0 
approximate formula in the time domain, that may be 
applied to arbitrary, non-homogeneous and anisotropic, 
albeit ohmic volume conductors: 


V()~ | a(t7)-Bs(@F)dv (10) 
B 


In this formula (a direct extension of Equation (1)) the 
effect of the non-uniform and transient polarization of 
the double layer of the active electrode on the lead field 


= = 1 Ese => i : 
u(t,r)= FIO?) is included introducing an 


explicit time dependence that must be calculated from 
the electric, geometric, and electrochemical properties 
of the system. 

But for many tissues, and in the case of myocardium in 
particular, the volume conductor is anisotropic and 
heterogeneous, and certainly bounded, so the study of 
the relation between the size and shape of the 
electrode’s head and the generated lead field u(t,7), 
including the location and geometric dimensions of the 
critical region for sensing bioelectric signals, should be 
extended to embrace this more general case. 
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Introduction of an auxiliary virtual space 
(alibi) 


For common types of heterogeneity and anisotropy, like 
the ones found in ventricular myocardium, a change of 
spatial variables from Euclidean to curvilinear 
coordinates may be found, allowing the reduction of 
several complex three- dimensional propagation 
problems of action potentials to propagation in the 
homogeneous and isotropic case, if an auxiliary 
“virtual space” with a smooth one to one 
correspondence with physical space, is introduced [14]. 

It is possible to apply this idea of an auxiliary space 
(alibi), to the description of the sensing sensitivity field 
(lead field) of an active electrode in a volume 
conductor. 

In the auxiliary space, Laplace’s equation for the 
potential and all the previous developments are valid, 
including the approximate analytical formulae relating 
the geometrical dimensions of the critical region with 
the distribution of electric potential along the axis of 
symmetry of the concave active electrode given in [1]. 
Going back to the original coordinates, the results are 
interpreted in physical space and the modifications that 
the lead field and the sensing sensitivities suffer due to 
the anisotropy and heterogeneity of the tissues are 
assessed. 

Let us briefly sketch the procedure by a typical 
example: the myocardium. (All the necessary 
background information that is needed from now on 
may be found in [15]). 

We assume that a small sensing electrode is localized in 
the ventricular cavity, touching the endocardium. 

The endocardial surface may be represented, 
approximately enough, as a plane. 

The myocardial fibres run almost parallel in layers, but 
the direction of the fibres in each layer varies smoothly 
from the endocardium to the epycardium. 


We introduce first a trihedralT, of unit vectors, 
hatoske with the origin and the unit vectors i, jo 


fixed in the endocardial plane, and ko orthogonal to this 


plane and pointing towards the epycardium. 
Now, taking as new origin each point of coordinate Z of 


the line defined by ko and the origin of T,, we 
introduce a new orthogonal trihedral T, of unit vectors 
@, =€ 5 ,€),€3 = ky, where é ; gives the direction of the 


fibres in the layer located at a depthz from the 
endocardium. 

Then any point in the myocardium may be given by the 
coordinate Zand its three coordinates X,,X,,X; 
relative to T,. 

The following connection exist between the coordinates 


x,y relative toT, and x,, xX, relative to T, : 


x = xsinla Z |+ycos|a Z (11a) 
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% =—xcos|a Z |+ysin|a Zz (11b) 


Here a(z) is the angle between the direction of the 


fibres in the layer, é - and the unit vector Te : 


IfGis the conductivity 
potential, then: 
V-(G-VV)=0 (12) 


Now, the conductivity tensor of the volume conductor is 


tensor andV the electric 


diagonal relative to T, in the corresponding layer of 
myocardial fibres. 

With enough approximation, its diagonal components 
G11, G5, G33 are functions of only one coordinate: 


Gi (x2), Gro(x2), Ga3(xs) 
As consequence, if V(x,,x,,x3) is the representation of 


the electric potential in local coordinates, equation 
(12) may be written thus: 
2.0 0 
AGG) Z—V1=0 (13) 
WON Ox, 


Nevertheless, we are still in physical space, naming the 
same points in another form (alias, using other 
coordinates). 

Now, let us introduce a virtual space (alibi) of 
coordinates €,,&,,&,, filled by a homogeneous and 
isotropic volume conductor of arbitrary (scalar) 
conductivity Go . 

The connection between corresponding points of each 
space is given by: 


XK 


G, 
fa d: k= 
sr I Gy (s) : 


The connection between the local representation of the 


1,2,3 (14) 


potential in physical space V(x,,x,,xX;) and the 
potential in the corresponding point of virtual space 
virtual space #(€,, €,, &) is defined by: 


WG E553) =V (4 (Gy (E5), 3(&)) 
(15) 
Furthermore, we introduce the functions: 
Siu) = Gy le (G)) =1,23 (16) 


The, from (14), (15) and (16), the left-hand member of 
(13) may be shown to be equal to: 


1 dg,. Ow 
G, Pa ( PG). Ga) 
oe 7 Bi AE, Of, 
So (13) is equivalent to: 
3. 92 F 
eae = tH) (18) 
tal OF; 2 Ly ix 
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ey 


O 

Here: Ly, yy & Se (19a) Le ga * ra (19b) 
( 
Og? de, 


If the spatial scales of the components of the gradient of 
the potential (the components of the electrical field) in 
virtual space L,, ;, are (as a rule) at least an order of 


magnitude less than the spatial scales of the 
components of the conductivity tensor in virtual space 


La ak , then the potential in the auxiliary space 
w(E,,6,3) will verify, approximately Laplace’s 
equation: 

ay 
k=l Ge, 7 


The approximation will be better the smaller are the 


Me 


(20) 


quotients of the spatial scales L,, ;,.and L eke * 


Conclusions 


(a) An experimental procedure for a parametric 
characterization of the lead field of sensing electrodes 
(concaves or not) from experimental data was 
described. 


Giving the good adjustment obtained by the application 
of Levenberg-Marquardt algorithm, and introducing the 
experimental values of 7, £,, and f,, then the 


components of the lead field of a polarizable electrode, 
in steady state and immersed in an electrolytic solution, 
could be predicted using formulae (5), at least near the 
axis of symmetry where the critical region of higher 
sensing sensitivity of concave electrodes is located. 
This kind of experiment in electrolytic tanks, with 
active electrodes of scaled sizes, and with suitable 
surface microstructures, allows an easy determination 
of the lead fields and the corresponding sensing 
sensitivities when the double layer is polarisable and 
polarizes non-uniformly, like in concave electrodes. 


(b) For unbounded, isotropic, and homogeneous volume 
conductors, the time dependence of the lead field of a 
polarizable active may be traced back to an explicit 
time dependence electrode of the coefficients of the 
multipolar expansion of the electric potential. 


The multipolar coefficients #, and £, can then be 


related, at its turn, with the size and shape of the 
sensing electrode’s head. 


Then, formula (10) can be used to relate the voltage 
signal in each lead with the distribution of biological 
electromotive forces that produce it, for measurement 
chains of input impedance at least an order of 
magnitude greater than the impedance of the lead 
system itself. 


An extension of formula (10) to the case in which this 
last assumption doesn’t apply would be of interest. 
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(c) The design problem of determining the geometric 
parameters of an electrode’s head such that it produces 
a critical sensing region for bioelectric signals, with 
suitable shape and size, as was posed and solved in [1] 
for homogeneous and isotropic volume conductors, can 
be extended to heterogeneous and anisotropic volume 
conductors introducing a virtual space that allows a 
reduction of the problem to the formerly solved one. 
After that the results can be translated to the 
heterogeneous and anisotropic case in physical space. 


(d) As the electric conductivity of the electrode’s 
material is at least several orders of magnitude greater 
than the conductivity of the tissues and the body fluids, 
the surface of the electrode may be safely assumed to 
be a surface of constant value of the electric potential 
field V produced when the reciprocal current is 
injected to the lead from outside. 


In principle any surface V =constant in the virtual 
space could be a candidate to a surface of a suitable 
electrode, and for each one of this, we have a 
corresponding surface in physical space. 


However, this cannot be exactly this way if the surface 
of equal potential is unbounded in physical space. 


But if our interest is focussed into the abovementioned 
critical region located near the axis of symmetry, it is 
possible to construct a suitable approximation to the 
shape and size of the electrode to produce the desired 
lead field and thus the desired distribution of local 
sensitivities within established tolerance bounds, as 
shown in [1]. 


(e) Now, it seems that enough groundwork has been 
done to enable an efficient application of numerical 
codes for elliptic partial differential equations to design 
problems posed by the sensing of bioelectric signals by 
concave electrodes. 


Nevertheless, the digital simulation of the non-uniform 
polarization process coupled with the field equations 
poses several problems. 


This, as well as specific examples of design of concave 
electrodes for sensing bioelectric signal will be the 
purpose of the last paper in this series. 
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